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1.  INTRODUCTION 


A  variety  of  problems  of  physical  interest  involving  the  deformation  of  elastic- 
plastic  solids  may  require  the  admission  of  discontinuities  in  the  gradients  of 
stresses  and  velocities  (weak  discontinuities)  or  in  these  quantities  themselves 
(strong  discontinuities).  Such  discontinuities  may  occur  within  regions  which 
are  currently  deforming  plastically  or  at  elastic-plastic  boundaries.  These  possi¬ 
bilities  have  received  wide  attention  for  rigid  -  perfectly  plastic  solids  in  plane 
strain  (Hill  [3])  and  in  generalized  plane  stress  (Hill  [l])  in  the  presence  of 
either  the  isotropic  Huber-von  Mises  or  Tresca  yield  conditions  in  the  plastic 
range.  It  is  well-known  that  for  such  solids,  strong  discontinuities  in  stress  and 
velocity  cannot  be  simultaneously  present,  and  that  velocity  jumps  occur  across 
characteristic  surfaces.  It  has  been  noted  by  Hill  [l]  that  when  a  rigid-plastic 
generalized  plane-stress  theory  is  employed  m  the  study  of  the  extension  of  thin 
plates,  two  types  of  strong  discontinuities  must  be  considered.  These  arise 
because  of  the  neglect  of  elastic  deformation  and  the  averaged  nature  of  gen¬ 
eralized  plane  stress.  A  consideration  of  the  second  of  these  factors  has  led  to 
the  mathematical  idealization  of  the  experimentally  observed  phenomenon  of 
localized  necking  in  thin  sheets  (Nadai  [8]). 

In  a  recent  paper  [6],  Drugan  and  Rice  investigated  strong  discontinuities 
across  quasistatically  propagating  surfaces  in  elastic-plastic  solids  under  gen¬ 
eral  three-dimensional  conditions  when  all  displacement  components  are 
assumed  to  be  continuous.  One  important  conclusion  of  the  work  reported  in 
[6]  is  that  all  stress  components  are  always  continuous,  a  result  that  follows 
from  certain  material  stability  postulates. 

Pan  [7]  has  also  discussed  quasistatically  moving  strong  discontinuities  for 
elastic-perfectly  plastic  Huber-von  Mises  materials  under  generalized  plane 
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stress.  He  assumes  that  a  strong  discontinuity  can  be  replaced  by  a  transition 
layer  of  elastic  material  in  which  all  stress  components  are  assumed  to  vary 
continuously.  He  subsequently  demonstrates  full  continuity  of  stresses  across 
propagating  surfaces,  by  using  the  specific  nature  of  the  Huber-von  Mises  locus 
and  arriving  at  a  contradiction. 

In  the  present  work  quasistatic  discontinuities  are  reexamined  for  the  more 
general  case  of  an  anisotropic  hardening  solid  using  an  integral  form  of  the 
maximum  plastic  work  inequality  and  the  usual  assumptions  in  the  theory  of 
generalized  plane  stress  (Section  2).  It  is  demonstrated  in  Section  5  that  the 
use  of  the  maximum  plastic  work  inequality  leads  to  full  stress  continuity  for  a 
broad  class  of  solids  which  includes  some  hardening  materials  and  anisotropic 
behavior.  Pan's  assumptions  and  the  limitations  of  his  approach  are  discussed 
in  Section  5.  A  complete  analysis  of  all  possible  velocity  jumps  including  sliding 
discontinuity  and  localized  necks  is  carried  out  in  Section  6  with  some  general¬ 
ity  in  constitutive  behavior.  These  results  are  then  specialized  for  the  elastic- 
perfectly  plastic  Huber-von  Mises  solid.  A  simplified  expression  is  obtained  for 
positive  plastic  work  accumulation  due  to  the  passage  of  a  discontinuity  surface. 

One  important  application  of  the  results  of  this  work  is  to  quasistatic  crack 
growth  problems  in  elastic-plastic  solids  under  conditions  of  generalized  plane 
stress  (Rice  [9]).  The  stress  and  velocity  fields  near  the  tip  of  a  growing  crack 
may  be  constructed  by  assembling  plastic  and  elastic  unloading  sectors  and  by 
satisfying  proper  matching  conditions,  boundary  conditions  and  material  stabil¬ 
ity  postulates.  Some  restrictions  arising  from  the  jump  conditions  on  this 
assembly  are  pointed  out  in  Section  7. 
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2.  THE  GENERALIZED  PLANE-STRESS  PROBLEM 

Consider  an  elastic-plastic  body  occupying  an  open  cylindrical  region  R  of 
height  2h  (see  Figure  l).  , 


Figure  1  Elastic-plastic  body  with  discontinuity  surface  E. 


Let  the  boundary  9R  of  the  above  region  be  composed  of  two  traction  free 
planar  surfaces  Sj  and  Sg  and  a  lateral  surface  L. 

Consider  further  a  fixed  orthonormal  coordinate  system  £o',e  i,e  2,e  3!  such 
that  e  3  is  parallel  to  the  generators  of  R. 

Generalized  plane  stress  conditions  require  that  the  height  of  the  cylinder 
(also  referred  later  as  the  thickness  of  the  cylindrical  plate)  be  small  as  com¬ 
pared  with  any  other  dimension  of  the  cylinder,  and  that  the  prescribed  trac¬ 
tions^  be  such  that: 


t  =  0  or  a31  =  0  on  Si  and  S2 


and 


t3  =  0.  ta  =  ta(x!,X2)  on  L, 


(2.1) 


Here  ay  are  the  components  of  the  symmetric  Cauchy  stress  tensor,  Greek  sub¬ 
script  have  the  range  1,  2  while  Latin  subscipts  take  the  values  1,  2,  and  3.  (This 
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convention  will  be  adopted  through  the  following  development.) 

In  what  follows  field  quantities  such  as  a  ,  e  ,  u  ,  and  v  will  represent  thick¬ 
ness  averages  of  the  stress  and  strain  tensors  and  the  displacement  and  velocity 
vectors  respectively.  It  is  also  assumed  that, 

ff3j  =  0  on  R.  (2.2) 

The  above  assumptions  result  in  solutions  of  the  generalized  plane  stress 
problem  which  in  general  will  not  satisfy  the  exact  three  dimensional  field  equa¬ 
tions  as  discussed  in  detail  by  Timoshenko  [5],  and  Hill  [3].  This  is  because 
some  of  the  compatibility  equations  are  not  generally  satisfied  and  errors  are 
involved  in  using  the  averaged  quantities  in  the  constitutive  law  and  the  yield 
condition.  However  if  the  plate  thickness  is  sufficiently  small,  the  generalized 
plane  stress  solution  is  expected  to  provide  an  accurate  approximation. 

Let  E  be  a  planar  surface,  parallel  to  the  X2-X3  plane  dividing  the  region  R  in 
two  open  subregions  R+  and  R~  such  that 

R  =  R+  u  R"  u  E 

We  will  define  the  normal  n  (x°)  to  E  at  a  point  x°  e  E  as  the  outward  normal  of 
the  closed  subregion  R"(R±  =  R*  u  E)  at  the  same  point  x° 

In  what  follows  the  surface  E  will  be  viewed  as  a  potential  surface  of  strong 
discontinuities  (discontinuities  in  stresses  and  strains)  and  will  be  allowed  to 
translate  quasistatically  with  a  normal  velocity  Vn  . 

Since  the  approximate  theory  of  generalized  plane  stress  treats  the  thickness 
of  the  plate  as  vanishingly  small,  Hill  [1]  points  out  that  every  quantity  whose 
gradient  is  of  order  (l/h)  in  a  zone  of  breadth  comparable  to  h  should,  be 
modeled  as  a  discontinuity.  Thus  the  experimentally  observed  formation  of 
necks  (Nadai  [8])  in  thin  plates  subjected  to  tension  (rapid  variation  of  thick- 
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ness  of  the  plate  in  narrow  zones)  would  be  modeled  as  discontinuities  in  the 
out-of-plane  displacement  component  Ug. 


The  jump  in  a  field  quantity  g  (x  ),  across  the  surface  E,  will  be  denoted  by: 


[&  ]  s  (x  )  ^ere 

g±  (x)  =  11m  g(x°  ±  sn  (x°))  x°  e  S  and  e>  0 


(2.3) 


3.  SMOOTHNESS  CONSIDERATIONS 

All  field  quantities  will  be  referred  to,  with  respect  to  an  orthonormal  frame 
$0,e  j,e  2,e  3]  translating  with  the  surface  E  and  such  that  OeS,  e3  =  ea  and 
e  i=n  ;  see  Figure  2. 


R 


%z 


R 


—2 

Figure  2.  Local  coordinate  system  translating  with  the  surface. 


In-plane  displacement  components  u„  are  required  to  have  the  following 
smoothness  properties: 

Ug  e  C(R)  and 

UgEC^R-E)  <'3'1^ 


with  the  understanding  that 


3Ug_ 

9x» 


need  not  be  continuous  across  E.  Then 


according  to  the  Hadamard  compatibility  relations  [2]  for  jumps  in  the 


derivatives  of  a  continuous  function, 
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=  \a  n^  on  E  (3.2) 

where  Xa  are  arbitrary  functions  of  position  on  E.  The  out-of-plane  displace¬ 
ment  component  u3  will  in  general  be  allowed  to  suffer  a  jump  across  E,  as  dis¬ 
cussed  in  Section  2.  Thus: 

Ug  e  C1(R+)  u  Cl{R~)  (3.3) 

with  the  understanding  that  on  E,  u3  and  its  gradient  need  not  be  defined.  On 
the  other  hand,  [u3],  the  jump  in  the  limiting  value  of  u3  from  R"  to  R+  will  be 
assumed  to  be  a  continuous  and  continuously  differentiable  function  of  position 
on  E. 


dn^ 

fop 


It  is  now  possible  to  extend  the  Hadamard  compatibility  relations  (3.2)  for 
the  treatment  of  jumps  in  the  derivatives  of  discontinuous  functions.  This 
extension  was  first  discussed  by  Thomas  [4],  The  following  simpler  version  was 
later  provided  by  Hill  [2], 


du3 

9Xj 


=  \3  nj  + 


on  E 


(3.4) 


where  X3  is  an  arbitrary  function  of  position  on  E  and  <p  is  an  arbitrary  continu¬ 
ous  function  together  with  its  gradient  on  E  and  in  one  neighborhood,  say  R+, 
with  the  additional  restriction  that 


9  -  [u3]  on  E. 

One  choice  of  <p  in  R'1'  would  be  to  consider  <p  continued  analytically  along  the 
normals.  Any  other  choice  would  merely  change  Ag  which  is  given  by 

Ag  =  [  ¥  u3  •  n  ]  -V  cp  ■  n 

Relations  (3.2)  and  (3.4)  allow  definition  of  jumps  in  the  strains  across  E  con- 
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sistent  with  the  assumptions  of  the  approximate  theory  of  generalized  plane 
stress. 

Within  the  contexts  of  a  small  strain  formulation, 

SH  =  jKi  +  uij)  (3-5) 

and  the  jumps  in  the  in-plane  strain  component  can  be  expressed  by  (3.2)  as: 

M  +  V1*)  orL  2  (3-6) 

On  the  other  hand,  the  jump  in  the  out  of  plane  strain  component  s33  can  be 
expressed  by  (3.2)  and  (3.4)  as: 

[£3s]  =^3n3  +  on  ^  .  (3.7) 

where  <p  e  (^(R*)  and  $p=[u3]  on  S. 

4.  MATERIAL  IDEALIZATION 

Within  the  context  of  the  small-strain  flow  theory  of  plasticity,  the  total 
strain  rate  tensor  can  be  decomposed  into  elastic  and  plastic  parts: 

£=£re  +  sp  onR  (4.1) 

r\j  ~ 

where  the  dot  denotes  differentiation  with  respect  to  time.  The  elastic  strain 

rate  tensor  is  related  to  the  stress  rate  tensor  <x  through  a  constant,  posi¬ 

tive  definite  four-tensor  5  (the  inverse  of  the  elasticity  tensor  C  ).  5  is 
assumed  to  possess  the  usual  major  and  minor  symmetries.  For  an  anisotropic 
elastic-plastic  solid,  s  e  is  given  by: 

se  =  S  a  on  R  (4.2) 

rsi  f-fci 

The  condition  for  the  existence  of  non-zero  plastic  strain  rates  k  p,  called  the 
yield  condition,  is  formulated  on  the  basis  of  a  scalar  valued  function  of  the 
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stresses  and  the  total  accumulated  equivalent  plastic  strain  (e11  )  denoted 
by  f(d  ,1?)  and  on  an  appropriate  definition  of  loading,  unloading  or  neutral 
loading.  The  condition  f(a  ,  ep)=0  defines  a  hypersurface  in  6-dimensional  stress 
space,  called  the  yield  surface.  In  this  work,  attention  will  be  focussed  on  the 
particular  class  of  materials  obeying  Drucker’s  stability  postulate.  A  particular 
form  of  Drucker’s  postulate,  known  as  the  maximum  plastic  work  inequality  can 
be  expressed  as: 

(a  -  a  *)  •  i  p  ^  0  (4.3) 

Vf (a  ,£P)=0,  and  f(<j  ’,1?)^  0.  The  two  most  important  implications  of  the  above 
postulate  are  the  following: 

a.  The  hypersurface  f(a  ,£P)=0  is  convex. 

b.  The  plastic  strain  rate  e  p  is  normal  to  the  yield  surface,  and  the  flow  rule 
takes  the  following  form: 

sp  =  \P  (4.4) 

where  A>  0  and  P  =  V  „  f.  k  and  P  are  scalar  valued  and  symmetric  tensor 
valued  functions  of  a  respectively.  In  the  following  section,  an  integral  form  of 
(4.3)  will  be  used  in  conjunction  with  equations  (4.1)  and  (4.2)  as  well  as  the 
compatibility  conditions  for  the  jumps  in  total  strains  (3.6),  (3.7)  to  define  the 
jumps  in  the  stresses  and  the  plastic  strains  produced  during  the  passage  of  a 
discontinuity  2  through  a  material  point. 

5.  STRESS  CONTINUITY  ACROSS  THE  PROPAGATING  SURFACE 

In  this  section  it  will  be  demonstrated  that  all  stress  components  sire  con¬ 
tinuous  across  the  surface  2  propagating  quasistatically  through  the  thin  plate. 
It  will  be  shown  that  this  is  true  even  if  the  out  of  plane  displacement  U3  suffers 
a  discontinuity  across  2.  The  following  proof  is  based  on  the  maximum  plastic 
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work  inequality  and  the  positive  definiteness  of  S  ,  It  is  an  extension  of  the 
proof  given  by  Drugan  and  Rice  [6]  for  the  general  three-dimensional  case.  In 
the  present  analysis,  only  the  in-plane  displacement  components  ua  are 
assumed  continuous,  and  the  proof  is  adapted  to  suit  the  assumptions  of  the 
theory  of  generalized  plane  stress.  Also,  unlike  the  discussion  by  Pan  [7]  and 
consistent  with  the  assumptions  of  generalized  plane  stress,  [l]  our  discussion 
treats  necks  as  jumps  and  not  as  narrow  transition  layers. 

If  inertia  terms  are  neglected,  the  balance  of  linear  momentum  requires  that 
across  the  quasistatically  moving  surface  E  the  traction  should  be  continuous. 
Thus 


[  ta]  =  |>af)  ]  =  0  on  2. 

With  respect  to  the  local  orthonormal  coordinate  frame  $0,e  i,e  2,e  3J  moving 
with  2,  n^djj  and  the  above  conditions  become: 

[crla]  =  0  on  2  (5.1) 


Equations  (2.2)  and  (5.1)  imply  that  the  only  stress  component  that  can  suffer  a 
non-trivial  jump  is  oZZ-  The  plastic  work  Wp  accumulated  discontinuously  at  a 
material  point  due  to  the  passage  of  the  surface  2  is  given  by: 


Wp  =  f  a  d  e  ’ 

/p+  ~  ~ 


(5.2) 


It  should  be  observed  here  that  some  error  is  involved  in  using  the  averaged 
stress  and  strain  quantities  of  generalized  plane  stress  in  the  above  integral. 
The  above  integral  is  evaluated  according  to  the  assumptions  of  Section  2.  On 
applying  equation  (2.2),  the  plastic  work  accumulation  in  (5.2)  reduces  to, 
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(5.3) 


Using  (5.1)  the  above  becomes: 

c22 

Wp  =  -CTu[ef i] -2ff12[£f2]  +  ^  a  22ds$2  (5.4) 

c23 

Also,  by  using  the  fact  that  n^  =  <5^,  (3.6)  implies  that: 

[e22]  =  0  on  E  or  [e|2]  =  -  [e22]  onE  (5.5) 

By  setting  ds  p  =  d  £  -  ds  8,  using  the  continuity  of  s22  across  E  (5.5)  and 
integrating  by  parts,  (5.4)  becomes 

C22 

Wp  =  -a n[ePi]  -  2 C712[ef2]  -  /  cr22d£|2.  (5.6) 

The  integral  in  (5.6)  can  now  be  evaluated  by  using  equations  (2.2),  (5.1),  and 
the  constitutive  law,  to  give: 

e|2 

f  ^aadsla  =  -  g-S2222(CT22  +  a22)[cr22]  (5.7) 

e22 

In  addition,  from  (2.2),  (5.1),  and  (5.5), 

[epa]  =  ~  [^la]  =  -  Saaaa  [ffaa]  011  £  (5.B) 

Thus  equations  (5.6)  and  (5.7)  give: 

WP  =  -  ff  ii[ePi]  -  2(7 !2[eP2]  -  2^022  +  a 2z)  [£aa] 
or, 

Wp  =  -|-(aif  +  a1r)[^].  (5.9) 

The  above  equation  is  of  the  same  form  as  the  corresponding  result  obtained  by 
Drugan  and  Rice  [6]  for  the  three-dimensional  case.  It  should,  however,  be 
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borne  in  mind  that  unlike  the  equivalent  expression  in  [6],  (5.9)  was  obtained  by 
requiring  that  only  the  in  plane  displacement  components  are  continuous  and 
by  using  the  assumptions  of  generalized  plane  stress.  It  should  also  be  observed 
that  the  restrictions  imposed  on  the  path  in  stress  space  in  the  evaluation  of 
the  integral  in  (5.6)  are  the  plane  stress  conditions  and  the  continuity  of  trac¬ 
tions  across  E.  This  effectively  implies  a  straight  line  path  in  stress  space  from 
O’ 22  tO  a 22- 

The  integral  form  of  the  plastic  work  inequality  (4.3)  can,  now  be  used  by  set¬ 
ting  a  *=ct  +  where  f(<7  '*',£P)  <  0.  Thus,  a  +  is  constrained  to  remain  always  at  or 
inside  the  yield  surface  during  passage  of  E.  Thus  by  using  (2.2),  (5.1),  and  (5.9): 

(cry  -  (7jj')dEjj>  =  -  +  &&)[*&]  +  ^  0 

which  by  (5.8)  gives 

i|“  [<7zz]  [£22]  -  0  or  ([022])®  S2222^  0  (5.10) 

(5.10)  now  requires  that  [022]  -  0  since  S2222  >  0. 

Remarks 

The  following  remarks  are  relevant: 

1.  Under  generalized  plane  stress  conditions,  all  stress  components  are  con¬ 
tinuous  across  the  slowly  propagating  surface  E  even  if  the  out  of  plane  dis¬ 
placement  U3  suffers  a  discontinuity.  Drugan  and  Rice  [6]  establishes  full 
stress  continuity  under  general  three-dimensional  conditions,  with  all  dis¬ 
placement  components  assumed  continuous. 

2.  The  present  discussion  applies  to  general  anisotropic  elastic-plastic  harden¬ 
ing  solids  obeying  a  flow  rule  of  the  associated  type.  The  proof  of  full  stress 
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continuity  is  based  on  an  integral  form  of  the  maximum  plastic  work  ine¬ 
quality  and  the  positive  definiteness  of  the  elastic  potential.  The  initial  yield 
surface  can  be  anisotropic  but  should  be  convex.  All  types  of  hardening  in 
which  the  subsequent  yield  locus  fully  encloses  all  previous  yield  loci  can  be 
included. 

3.  An  earlier  discussion  by  Pan  [7]  is  limited  to  elastic-  ideally  plastic  solids  of 
a  Huber-von  Mises  type  under  generalized  plane  stress  conditions.  His  argu¬ 
ment,  which  does  not  make  use  of  the  maximum  plastic  work  inequality,  fol¬ 
lows  from  Hill’s  [3]  statement  that  the  stress  state  from  cr22  to  a can  be 
bridged  only  by  a  succession  of  elastic  states.  This  assumes  a  smooth  varia¬ 
tion  of  stresses  in  a  'transition  layer."  Such  an  assumption  is  questionable 
for  generalized  plane  stress  since,  as  pointed  out  in  Section  2,  any  field 

quantity  whose  gradient  is  O(^j-)  in  a  zone  of  breadth  comparable  to  h 

should  be  modelled  as  a  discontinuity.  Even  if  this  assumption  is  accepted, 
Pan’s  argument  clearly  does  not  apply  to  arbitrary  yield  surfaces  or  general 
hardening  solids.  For  instance,  in  elastic-perfectly  plastic  solids  character¬ 
ized  by  a  Tresca  yield  condition  when  the  neck  (discontinuity  in  u3)  coincides 
with  a  principal  stress  direction  and  an  =  i  a0,  the  stress  component  a2 z 
can  have  any  value  between  0  and  ±  a0  and  still  lie  on  the  yield  surface  (Hill 
[3]).  Hence  A  in  equation  (4.4)  is  not  necessarily  zero  in  the  transition  from 
ct  +  to  a  ”  (Hill  [l])  and  the  argument  fails.  Also  for  any  type  of  hardening 
solid,  the  consistency  condition  requires  the  stress  state  to  lie  on  the  yield 
surface  during  the  process  from  <r  *  to  ct  “  and  no  elastic  unloading  is  possi¬ 


ble. 


-  13  - 


6.  DISCONTINUITIES  IN  STRAINS  AND  VELOCITIES 

In  this  section,  the  earlier  result  pertaining  to  continuity  of  stresses  across  2 
will  be  used  to  provide  restrictions  on  the  nature  of  admissible  jumps  in  strains 
and  material  particle  velocities  across  E  for  a  general  anisotropic  hardening 
solid.  Attention  will  then  be  turned  to  plastically  incompressible,  generally 
anisotropic,  elastic-perfectly  plastic  solids  with  smooth  but  otherwise  arbitrary 
yield  surfaces.  Specialized  results  will  be  given  for  Huber-von  Mises  solids  at  the 
end  of  the  discussion. 


General  Considerations: 


The  jumps  in  the  in-plane  velocity  component  va  are  given  ([2],  [6])  by: 

M  =  (6.1) 

where  Vn  is  the  normal  velocity  of  E.  Making  use  of  (3  2)  and  (3.6)  the  velocity 
jumps  may  be  expressed  as: 


[Vi]  =  -VOn] 

W  =  -2V[eia]  °n  2 


(6.2) 


Full  stress  continuity  and  (4.2)  require  the  elastic  part  of  the  strains  to  be  con¬ 
tinuous  across  E 


[ejj]  =  0  on  E  (6.3) 

The  above,  and  equation  (5.8),  therefore  imply 

[efg]  =  0  on  E  (6.4) 

As  a  result,  the  expression  for  the  positive  plastic  work  accumulation  in  (5.9) 
becomes 


W1*  -  -  o  n[«u]  —  lafefa]  —  0 


(6.5) 
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and  the  jumps  in  the  velocity  components  V!  and  v3  are  given  by: 


[▼i]  =  -V[eFi] 

[v2]  =  -  2V  [ef2] 


(6.6) 


The  plastic  work  Wp  can  now  be  expressed  in  terms  of  velocity  jumps  as  follows: 


Wp  =  ^K<Tn[v,]  +  a is[v2])  >  0 


(6.7) 


No  specific  restrictions  on  the  constitutive  model  other  than  the  general 
assumptions  made  in  Section  4  have  been  imposed  in  the  derivation  of  equa¬ 
tions  (6.l)-(6.7). 


For  the  specific  class  of  plastically  incompressible  solids: 

O33]  =  “  On]  ~  0a2]  on  E  (6.8) 

which  by  use  of  equation  (6,4)  simplifies  to: 

[fi3a]  =  ”  [£n]  on  2  (6.9) 

Equation  (6  9)  serves  to  determine  the  jump  in  the  out  of  plane  plastic  strain 
component  in  terms  of  the  jump  in  the  in-plane  plastic  strain  component  eft 
for  plastically  incompressible  solids. 

If  the  displacement  component  u3  happens  to  be  continuous  across  E  as  in 
[6],  then  £33  and  hence  would  also  be  continuous.  Equations  (6,8)  and  (6.6) 
will  then  imply  that  and  vx  should  also  be  continuous  across  E.  Thus  it  fol¬ 
lows  that  for  a  plastically  incompressible  solid,  when  the  surface  E  does  not 
coincide  with  a  neck  (jump  in  u3)  only  a  sliding  velocity  discontinuity  (jump  in 
v2)  is  permissible. 


Elastic’ Perfectly  Hastic  Solid 


For  such  solids,  the  yield  surface  is  represented  by 
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f(<7  )  =  0  on  R  (6.10) 

where  f(a  )  depends  symmetrically  on  a  and  a  T.  It  will  also  be  assumed  here 
that  the  yield  surface  is  smooth  (has  a  continuous  normal) 

Under  such  circumstances  the  flow  rule  takes  the  following  form: 

ep  =  \P  onR  (6.11) 

rsj  t sj 

where  A>  0  is  an  undetermined  scalar  function  of  position,  and 

P  (a)  =  V,  f(ff  )  onR  (6.12) 

is  a  symmetric  tensor-valued  function  of  a  .  Under  conditions  of  generalized 
plane  stress,  equations  (6.10)  and  (6.12)  should  be  used  in  conjunction  with  the 
constraint  (2.2).  Inside  regions  which  are  currently  deforming  plastically,  it  can 
be  shown  from  the  two  in-plane  equilibrium  equations,  the  yield  condition  and 
the  plane  stress  assumption  (2.2)  that  along  stress  characteristic  directions  the 
direct  components  of  Pa/}  should  vanish  (Hill  [3]). 

It  is  also  clear  that,  P^  should  be  continuous  across  £  from  assumed  smooth¬ 
ness  of  the  yield  surface  and  the  requirement  of  full  stress  continuity.  Then, 
from  (6.11)  the  jumps  in  the  plastic  strain  component  sp  becomes: 

[Cl  =  -V  pij  onS  (613) 

X‘ 

where  rj  =  J  dA>  0  is  an  undetermined  scalar  function  of  position  on  2. 

\+ 

Since  [£22]  =  0  across  2,  equation  (6.13)  implies  that  either  77= 0  or  P22=Q  or 
both.  If  77=0,  (6.13)  requires  all  strain  components  to  be  continuous.  Thus  the 
necessary  condition  for  non-trivial  jumps  in  strains  to  exist  across  2  is  that  P22 
should  vanish  on  2.  In  other  words,  2  should  coincide  with  a  stress  characteris¬ 
tic  direction  of  its  plastic  side 
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This  condition  is  less  restrictive  than  the  necessary  condition  for  non-trivial 
jumps  in  the  plastic  strain  components  derived  by  Drugan  and  Rice  [6]  when  all 
displacements  were  continuous  across  E.  The  corresponding  necessary  condi¬ 
tion  derived  in  [6]  states  that  P22=P33=P23=0  on  E. 

From  the  above,  the  following  important  observation  can  be  made: 


Consider  at  least  one  side  of  E  (which  coincides  with  a  neck,  say  R+)  to 

currently  be  deforming  plastically.  If  in  addition,  E  coincides  with  one  of  the 

stress  characteristic  directions,  say  direction  A  (see  Figure  3), 

Choracteristic 


k 


k 


Figure  3.  Moving  discontinuity  surface  and  characteristic  directions  of  the 
plastically  deforming  side. 


then  the  velocity  component  along  the  other  characteristic  direction  B,  is  con¬ 
tinuous  across  E.  Thus 


[Vb]  =  0  onE.  (6.14) 

The  above  follows  by  first  observing  that  since  E  coincides  with  a  stress 
characteristic  direction,  P22  vanishes  on  E.  Also,  if  the  other  characteristic 

direction  makes  an  angle  6  (6*  ±  —■)  with  the  xx  axis,  then  by  the  fact  that 

Pqb=0  and  the  transformation  relation,  we  have 
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p 

tan  <5  =  -  opn  for  P12  *  0.  (6.15) 

12 

In  addition,  combining  (6.6)  and  (6.13)  and  noting  that  P^O  the  following  is 
true: 

OJ  =  JFriv2]  on  2.  (6.16) 

*P  13 

The  velocity  jump  [vg]  along  the  other  characteristic  direction  will  be  given  by: 
[vB]  =  cos  5([vj]  +  [v2]  tan (5),  6  *  ± 

which  vanishes  by  use  of  (6.15)  and  (6  16).  This  general  result  was  also  noted  by 
Pan  [7]  for  the  special  case  of  an  isotropic  Huber-von  Mises  solid  and  it  also 
holds  for  stationary  necks  in  a  rigid-plastic  solid  (Hill  [1]). 

If  in  addition  P12=0,  both  the  stress  characteristics  merge  along  E  (5=  ±  ^-) 

and  as  a  result  E  becomes  a  '‘parabolic  line."  Equations  (6.6)  and  (6.13)  then 
imply  that  if  P12=0, 

[v2]  =  0  on  E.  (6.17) 

Thus  when  E  coincides  with  a  "parabolic  line,"  the  tangential  velocity  is  continu¬ 
ous  and  only  the  normal  velocity  has  a  jump. 

When  E  coincides  with  a  neck  and  the  two  characteristic  directions  do  not 
merge  along  E  (see  Fig.  3)  then  the  accumulation  of  plastic  work  (6.7)  due  to  the 
passage  of  E  becomes: 

wp  =  1~  (  OuPu_^12Pl2.)  [vg]  ^  0  (6.18) 

Also  the  fact  that  sfi  >  0  imphes  that 


on  P;j  >  0 


(6.19) 


-  18  - 


By  (8.2)  and  P22=0  along  E,  (6.19)  becomes: 

^iiPn  +2ai2P13>  0  (6.20) 

Inequalities  (6.18)  and  (6.20)  result  to 

p'2*0  <621) 

when  the  two  characteristics  merge  along  S  (<5=±:  and  it  follows 

from  (6.7)  and  (6.17)  that 

^-[vJ^O.  (6.22) 

Isotropic  Huber- von  Mises  solids 

The  above  results  can  now  be  specialized  for  an  isotropic  elastic-perfectly 
plastic  solid  which  obeys  the  Huber-von  Mises  yield  condition.  For  such  a  solid, 
the  yield  condition  states 

f(£)=  |-S  S  -  t|  =  0  onR  (6.23) 

where  S  =  a  -  ^-tr  o  1  is  the  deviatoric  stress  tensor  and  r0  is  the  yield  stress 

hu  Hu  3  Hu  hu  O  J 

in  pure  shear.  For  such  a  solid 

P  (£  )  =  V  2.  f(£  )  =  S  on  R  (6.24) 

All  the  results  and  corresponding  remarks  from  (6.3)-(6.22)  hold  for  this  solid 
with  P  replaced  by  S  .  In  particular  (6.16)  takes  the  form  (Pan  [7]): 

[Vi]  =  Sr[Vzl  if  Sia  * 0  (6-25) 


wP  =  v  (  0 

v  0 12 


and  (6.18)  reduces  to 
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(6.26) 


Sarrunary  of  Results 

The  results  of  Section  6  can  now  be  summarized  as  follows: 

a.  For  a  general  anisotropic  hardening  solid  which  is  also  plastically 
incompressible  the  following  is  true:  When  the  propagating  surface  £  does 
not  coincide  with  a  neck  (full  displacement  continuity),  only  a  jump  in  the 
tangential  velocity  component  (sliding  discontinuity)  is  admissible. 

b.  If,  however,  the  solid  is  perfectly  plastic,  £  coincides  with  one  characteristic 
direction  (P22=0)-  In  addition,  full  displacement  continuity  together  with 
plastic  incompressibility  also  give  P33=Pu=0.  This  states  that  the  direction 
normal  to  £  is  also  a  characteristic  direction.  Unlike  plane  strain,  this 
occurs  under  plane  stress  condition  only  under  exceptional  circumstances 
(Hill  [3]).  In  particular,  for  Huber-von  Mises  solids  this  is  true  when  the  sur¬ 
face  coincides  with  a  plane  of  maximum  shear  stress  and  the  latter  is  equal 
in  magnitude  to  the  yield  stress  in  pure  shear. 

c.  For  a  general  anisotropic  elastic-perfectly  plastic  solid,  when  a  surface  coin¬ 
cides  with  a  neck  (discontinuity  in  u3)  both  tangential  and  normal  velocities 
have  jumps.  This  requires  that  the  neck  should  lie  along  one  characteristic 
direction.  Then  the  component  of  the  velocity  along  the  other  characteristic 
direction  (not  generally  perpendicular  to  £)  is  continuous  (see  equation 
(6  14)).  Thus  necks  cannot  form  if  the  plastically  deforming  side  of  the  sur¬ 
face  is  in  an  elliptic  state  of  stress. 

d.  For  an  elastic  perfectly  plastic  solid  if  in  addition  to  (c),  P^O,  both  the 
characteristics  merge  along  the  neck  and  this  results  in  a  parabolic  stress 
state.  Then,  the  tangential  velocity  is  continuous  and  only  the  normal  veto- 
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city  has  a  jump.  For  the  special  case  of  a  Huber-von  Mises  solid,  Pi2=S12-0 
and  the  characteristic  surface  coincides  with  a  principal  stress  direction. 

7.  REMARKS  AND  APPLICATIONS 

The  jump  conditions  discussed  here  have  some  relevance  to  the  stress  and 
strain  fields  near  the  tip  of  a  quasistatically  growing  crack  in  an  elastic-plastic 
solid  under  generalized  plane  stress  conditions.  For  instance,  in  the  elastic- 
perfectly  plastic  Huber-von  Mises  material  (Rice  [9])  a  "constant  stress"  (asymp¬ 
totic)  plastic  sector  cannot  occur  directly  behind  a  "centered  fan"  plastic  sector 
because  the  condition  for  positive  plastic  work  accumulation  (6.21)  will  be 
violated  at  the  interface.  This  renders  the  asymptotic  solution  for  the  plane 
stress  stationary  crack  by  Hutchinson  [10]  unacceptable  when  the  crack  begins 
to  grow.  From  the  preliminary  asymptotic  analysis  by  Rice  [9],  it  then  follows 
that  only  an  "elastic  unloading"  sector  can  occur  behind  the  centered  fan. 
Hutchinson's  stationary  crack  solution  also  has  a  jump  in  the  in-plane  stress 
component  between  two  constant  stress  sectors.  This  is  also  inadmissible  when 
the  crack  begins  to  propagate. 

No  solution  for  this  problem  which  satisfies  all  the  conditions  set  forth  in  the 
present  paper  has  yet  been  constructed.  An  open  question  that  arises,  for  which 
detailed  experimental  and  numerical  studies  may  provide  an  answer,  is  whether 
necking  occurs  near  the  growing  crack  tip.  Otherwise,  except  in  special  cir¬ 
cumstances  (like  a  fan  angle  of  90°),  no  strong  discontinuities  near  the  growing 
crack  tip  can  be  admitted.  In  view  of  the  fact  that  the  (fully  yielded)  stationary 
crack  tip  solution  [10]  has  a  strong  discontinuity,  one  wonders  whether  the  con¬ 
dition  of  full  continuity  in  both  stress  and  velocity  near  the  propagating  crack 
tip  may  be  too  restrictive  to  satisfy. 
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8.  CONCLUSION 

Strong  discontinuities  across  quasistatically  propagating  surfaces  in  aniso¬ 
tropic  elastic-plastic  solids  under  generalized  plane  stress  have  been  reexam¬ 
ined  allowing  for  some  generality  in  constitutive  response  and  taking  into 
account  the  phenomenon  of  necking.  Jumps  in  stresses  have  been  ruled  out  on 
the  basis  of  material  stability  postulates  and  a  previous  approach  (by  Pan  [7]) 
has  been  discussed.  It  has  been  noted  that  for  elastic-perfectly  plastic  solids, 
sliding  velocity  discontinuities  occur  under  restrictive  and  exceptional  condi¬ 
tions  (when  both  the  surface  and  its  normal  are  stress  characteristics)  for  gen¬ 
eralized  plane  stress  as  compared  to  plane  strain.  Necks  may  form  along 
(stress)  characteristic  directions  with  the  relative  velocity  vector  orthogonal  to 
the  other  family  of  characteristics. 
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